STABILITY OF THE GAUGE EQUIVALENT CLASSES IN 
STATIONARY INVERSE TRANSPORT 
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' Abstract. For anisotropic attenuating media, the albedo operator determines the 

^SJ , scattering and the attenuating coefficients up to a gauge transformation. We show 

jj^' that such a determination is stable. 

^ ■ 1. Introduction 

P^ '. This paper concerns the problem of recovering the absorption and scattering prop- 

^ I erties of a bounded, convex medium C M", n > 2 from the spatial-angular mea- 

surements of the density of particles at the boundary dQ. Provided that the particles 
interact with the medium but not with each other, the radiation transfer in the 
steady-state can be modeled by the transport equation 



(1) 



- 9 ■ Vn(x, 9) - a{x, 9)u{x, 9) + [ k{x, 9', 9)u{x, 9')d9' = 0, 



■ for X G and 9 G S*""^; see, e.g. [101 ES]- The function u{x, 9) represents the density 

CN . of particles at x traveling in the direction 9, a(x, 9) is the attenuation coefficient at 

p 

in 



X for particles moving in the direction of ^, and 9\ 9) is the scattering coefficient 
(or the collision kernel) which accounts for particles from an arbitrary direction 9' 
^ ■ which scatter in the direction of travel 9. Let T± denote the incoming and outgoing 

g '. "boundary" 

(2) T±:= {{x,9) edVtx S""-^ : ±9-n{x)>Q}, 

^ . n{x) being the outer unit normal at a boundary point x G dVt. The medium is probed 

with the given radiation 

(3) u\r_ = /-. 

The exiting radiation u\r^ is detected thus defining the albedo operator A that takes 
the incoming flux /_ to the outgoing flux M|r+, i-e. A[f-] := M|r+- 

In general, the boundary value problem ([1]) and ([3]) may not be uniquely solvable 
but it has a unique solution under some physically relevant subcritical conditions like 
(125]) 5 or (I5T]) . We note, however, that for sufficiently regular coefficients, the 

problem has unique solution for generic (a, k), see [301 EB]- 

One of the inverse boundary value problems in transport theory is to recover the 
attenuation coefficient a and the scattering kernel k from knowledge of the albedo 
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operator A. This problem has been solved under some restrictive assumptions (e.g. 
A; of a special type or independent of a variable) in[Il[2lEl[Sl[ISl[ia[2niEIlE2]. In 
three or higher dimensions, uniqueness and reconstruction results for general k and 
a = a{x) were established in [12]. The general approach is based on the study of the 
singularities of the fundamental solution of ([1]) (see also [9]), and the singularities of 
the Schwartz kernel of A. Stability estimates for k of special type were established in 
[271 [33]; cind recently, for general fc, in [7|. Uniqueness and reconstruction results in a 
Riemannian geometry setting, including recovery of a simple metric, were established 
in [21]. Similar results for the time-dependent model were established in [11], and 
in [T3] for the Riemannian case. In planar domains the work in [22] shows stable 
determination of the isotropic absorption and small scattering, and an extension to 
simple Riemannian geometry is given in [22]. Also in two dimensional domains we 
point out that the recovery of k is only known under smallness conditions which are 
more restrictive than what is needed to solve the direct problem; e.g. more restrictive 
than fl35l) or fl36|) below. On the other hand, in the time-dependent case, the extra 
variable allows us to treat the planar case without such restrictions, see [11]. We 
also mention here the recent works [8l[T6l[T7j, in which the coefficients are recovered 
from angularly averaged measurements rather than from the knowledge of the whole 
albedo operator A. For an exhaustive account on the inverse transport problem we 
refer to the review paper [B]. 

The above mentioned results concern media with directionally independent absorp- 
tion a = a{x), except for transport with variable speed when the attenuation may 
depend on \v\, a = a{x, \v\). 

The attenuation accounts not only for the absorption of particles, but also for the 
loss of particles due to the scattering. In the physical case in which k depends on 
two independent directions, the attenuation is inherently anisotropic a = a{x,6). 
However, in an anisotropic attenuating media, the unique determination of the co- 
efficients from boundary measurements no longer holds: In [2B] it is shown that the 
albedo operator determines the pairs of coefficients up to a gauge transformation; see 
(jlj) below. This non-uniqueness motivates the following definition. 

Definition 1.1. Two pairs of coefficients (a, k) and (a, k) are called gaMge equivalent 
if there exists a positive G L~(fi x S"""^) with 6 ■ S/^(t){x,e) G L°°(fi x 5""^) and 
= 1 for X G dVl such that 

(4) a(x, e) = a(x, e)-e- v. iog0(x, e), k{x, e\ e) = ^^Kx, e\ e). 

We denote the equivalence class of{a,k) by {a,k), and the equivalence relation itself 
by ~. 

The relation defined above is reflexive since (a, k) ~ (a, k) via = 1; it is symmetric 
since {a,k) ~ {a,k) via yields {a,k) ~ {a,k) via 1/0; and it is transitive since if 
(a, k) ~ (a, k) via and (a, k) ~ (a', k') via then (a, k) ~ (a', k') via 00. 

The main result in [28] is that, in dimensions n > 3, ^ = ^ if and only if (a, k) ~ 
(a, k), i.e., uniqueness up to gauge transformations. The uniqueness up to the gauge 
transformation extends naturally to refractive media and to dimension two, see |23j . 
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In this paper we study the question of stabihty of the determination of the gauge 
equivalent classes. Let (M, || ■ \\m) and (A^, || • be Banach spaces in which the 
attenuation and, respectively, the scattering kernel are considered, {a,k),{d,k) G 
M X N. The distance A between equivalence classes with respect to M x iV is given 
by the infimum of the distances between all possible pairs of representatives. More 
precisely, 

(5) A{{a,k), {d,k)) := inf _ _ max{||a' — a'||jv/, ||/c' — ^'Htv}- 

{a' ,k')e(a,k),{a' ,k')e{d,k) 

For n > 3 we work within the class of coefficients 

(6) (a, k) E X 5"-^) X L^{n x 5""^ 

For two dimensional domains {n = 2) both coefficients are assumed bounded: 

(7) (a, k) e L~(fi X S^) X L°^{Q x x S^). 
The following norms are used throughout 

||a||oo = ess sup,^^^g-^^Q^sn-i\a{x,9)\, 

||/i;||oo,i = ess sup(^^g,)gs^xS"-i / \k{x,6',e)\de, 
\\k\\oo =ess sup(^^e,_^,)gj^^5ix5i \k{x,9',9)\, 
\\k\\^ = [ [ [ \k{x,9',9)\dxd9'd9. 

Note that the gauge transformations preserve the class of coefficients in ([6]) or (I7j). 

In Section [2] we reduce the original inverse problem in Q to the inverse problem of 
transport in a larger (strictly convex) domain D Q, where the attenuation and 
scattering coefficients are extended by zero in Br \ Q. More precisely we show that 
the difference of two albedo operators realizes an isometry when transported from 
dfl to dBji. For simplicity, the larger domain is a ball but this is not essential. Let 
(a, k) be another pair of admissible coefficients for which the forward problem in Q is 
well posed and let A denote the corresponding albedo operators. Set a = d = and 
k = k = in Br \ Vl. Then the forward problems in Br are also well posed and let 
and denote the corresponding albedo operators respectively. The boundary 
data is considered on 

(8) Vl:= {{x,9) eOBrX S""-^ : ±9-n{x)>Q}, 

n{x) now being the outer unit normal at a boundary point x G OBr. Provided that 
the forward problem is well-posed in L^, 1 < p < oo, we show that 

(9) M ~ -4|U(LP{r_;dO;-LP(r+;dO) = ~ -4^ll£(LP(r«;d§«);LP(r«;dC«))- 

In (|9]) we used d^ = \n{x) ■ 9\dfi{x)d9, where d9 is the normalized measure on the 
sphere, dfi{x) is the induced Lebesgue measure on dfl and n{x) is the unit outer 
normal at some x G dQ. Similarly, d^^ = \n{x) ■ 9\dfi^{x)d9, where dfi^{x) is the 
induced Lebesgue measure on OBr. 
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Consequently, we may consider the data (albedo operators) given directly on the 
dB^ and drop R from their notation. The isometry ^ with p = 1 is used for domains 
in three or higher dimensions. For brevity let 

11^ -^11 := M - ^||£(Ll(r_;d5);Ll(r+;dO)- 

Let r±(x, 9) be the travel time it takes a particle at x G Br to reach the boundary 
dBji while moving in the direction of ±^ and define r(x, 9) = r_(x, 9)+t+{x, 9). Since 
we work with unit-speed velocities, t{x,9) < 2R. Moreover, since dist(fi, 9-8/?) > 0, 
we have 

(10) Cr:=M{t{x,9) : (x, ^) G H x > 0. 

Note that we could make cr = 1 at the expense of a sufficiently large radius R. 
For domains in three or higher dimensions, and for S, p > we consider the class 

(11) f/s^p := {(a, /c) as in dZD : ||a||oo < S, ||A;||oo,i < p}- 
The main result of stability of gauge equivalent classes is the following. 

Theorem 1.2. Let {a, k), {a, k) G f/s.p be such that the corresponding forward prob- 
lems are well posed. Then 

(12) A{{a,k),{~a,~k))<C\\A-Al 

where A is with respect to L°° x , and C is a constant depending only onT,, p, cr 
and R. More precisely, there exists a representative (a', k') G (a, k) such that 

(13) ||a' - alloc < CP- ^11, 

(14) \\k' -k\\i<C\\A- A\\, 
where 

(15) C = max{7ri?e^«^ (l + 2pe^«^) , e^^^/c^}. 

For the stability of the equivalence classes in two dimensional domains we need a 
more refined notion of distance between the albedo operators: Following [291 Proposi- 
tion 1] , the Schwartz kernel of a albedo operator A admits the singular decomposition: 

A(x' 9') 

(16) a = ^^^^ ^ ^ S{^'+r+(x',9')0'}{x)S{0>}{9) + f3{x, 9, x', 9'), 



where 



/ rT+{x',e') 

(17) A{x', 9') = exp f - y a{x' + t9' , 9')dt 

and 1^ X 9'\(3 G L~(r^ x Ff!); see also Section [S 

Let y4, (3 be the coefficients corresponding to the decomposition f|T6|) of another 
albedo operator A. We define 

(18) P||. = max{||A|U;|||^x^'|/3||oo}. 
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In Section [2] we show that the isometric transportation from dQ to dBji holds in the 
new operator norm, i.e., 

^19^ II '■'4. t-^. 1 1 ^ — 1 1 1^4. 1^4. II** 

In the two dimensional case, for S,p > 0, we consider the smaller class 

(20) li;,p := {(a, fc) as in (I7j) : ||a||oo < S, ||A;||oo < p}- 

Theorem 1.3. For any S > 0, there exists < p < 1 depending only on S, R and 
cji, such that the following holds: If (a, k), (a, k) G V^.p, then 

(21) A{{a,k),{a,k))<C\\A-A\U, 

where A is with respect to x L°° , and C is a constant depending only on S, p, cr 
and R. 

The proofs of Theorem 11.21 and 11.31 are based on the analysis of the singularities of 
the Schwartz kernel of A as in [12] and, respectively, [29]. We also use an extended 
version of an estimate in [7], see Section |4] below. 

One can formulate and prove similar results in the case where the velocity belongs 
to an open subspace of M"", i.e., the speed can change, as in [71 [12]. We restrict 
ourselves to the fixed speed case (16*1 = 1) for the sake of simplicity of the exposition. 

2. Isometric transportation of the albedo operators from dVl to dBn 

Recall that Vt is strictly convex with Vt C Br, and A is the albedo operator for the 
radiation transport in Vt. The coefficients (a, k) are extended by zero in Br \ Vt and 
let Ar be the albedo operator corresponding to the radiative transport in Br x 5*""^, 
which takes functions on to functions on F^; see ([H])- Recall that t±{x,9) is the 
travel time it takes a particle at x G Br to reach the boundary SBr while moving in 
the direction of ±6'. 

We consider next the set of "projections" of T± onto F^ defined by 

(22) := {(x ± r±(x, 9)6, 9 : {x, 9) G F±} C F^ 
and the transportation maps 

(23) [rf]{x',9') := f{x'-T4x',9')9',9'), \/{x',9') G F_ 

(24) [ff]ix + T^{x,9)9,9):=f{x,9), V(x,e)GF+. 

T takes maps defined on to maps on F_, whereas T takes maps defined on F+ to 
maps on F;^. 

Proposition 2.1. For any 1 < p < oo, the maps 

(25) T:L^{r^-dO-^LP{r.;dO, 

(26) f:L^iT+;dO-L^iT^;dO 
are isomorphisms. 
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Proof. From the definitions fl23|) and fl2^ we have 

(27) ess suppii/ = ess supr_[T/] and ess suppii[T/] = ess supr_/. 
This proves the isometry for p = oo. For 1 < j> < oo we have the identities 

(28) / \f{x,e)\pdax,e)= [ \[Tf]{x',e)\pdax',e), 



(29) / \rf{x,e)\pdax,e)= \f{x',e)\pdax',e) 

obtained by the change of variables x = x' — t_(x', 6) and x = x' + t_|_(x', 6) respec- 
tively. □ 

Proposition 2.2. Let (a, k) be an admissible pair for the transport in Q such that 
the forward problem is well-posed in WiQ, x S*""^), 1 <p < oo, n > 2, and let A be 
the corresponding albedo operator. Extend the coefficients by zero in \ VL. Then 
the forward problem in U'{Bfi x S'^~^), 1 <p < oo, is also well posed and let be 
the corresponding albedo operator. For any f G L^(r:^) we have 



(30) A''[f]{x,d) 



f{x-T4x,e)e,e), if {x,e) er^\rl 
fATfix,e), if ix,e)eT^. 



Proof. The relation ( l30l) is a direct consequence of the fact that in \ Q, where 
the coefficients vanish, the solution u{x, 9) of the transport equation is constant in x 
along the lines in the direction of 9. □ 

Proposition 2.3. Let (a, k) and (a, k) be admissible pairs for the transport in Vt such 
that the forward problem is well-posed in L^{Q x S"'~^), n > 2, 1 < p < oo, and let 
A and A be the corresponding albedo operators. Extend the coefficients by zero in 
Br \ VL and let A^, respectively A^, be their corresponding albedo operators for the 
transport through B^. Then ^ and 0191) hold. 

Proof. Following (l30l) 

(31) {A"" - A'')[f]{x,9) = 1 ' ^ I , ; + \ +, 

\ f[A-A]Tf{x,9), zf {x,9)eTl 

For / G LP(r^), 

||(^^-i^)/|L.(r.) = ||t(^-^)T/||^,(f.) = ||(^-^)T/|U.(r,) 

< \\A-A\\ ■ ||T/|U.(r_) = \\A-A\\ ■ ||/|L.(f«) 

< M-^M|/|L.(r«). 

Hence ||^« - ^^|| < 
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To prove the converse inequality, let /o be the projection of / on r5, i.e. 



fo-.-- 

Then 



0, on \ 



/, on r 



R 



IIM-^)T/|U.(r,) = ||TM-^)T/|L,(f.) = ||M^-^^)/|L.(rH) 

= ||M^-^^)/o||L.(r?) < M^-^^ll ■ ||/o|L.(r«) 

= 11^ - ^ II ■ ||/|pfl||^p(pfl) = 11^ -A II ■ ||T/||z,p(r_). 

Since T is onto LP(r_), the inequality above yields ||^ — ^|| < ||^-^ — A^\\. 

Next, we prove the isometry in the || ■ ||*-norm. Let a,a,a^,a^ be the Schwartz 
kernels associated with the albedo operators A, A, A^, A^, and let A, A, A^, A^ and 
(3, (3, (3^, (3^ be the coefficients from the corresponding singular decomposition as in 

(USD. 

On the one hand, for / G C°°{T^) arbitrary, and (x, 9) G r+ we have: 
[fiA - A)rf]ix + r+(x, 9)9, 9) = [{A - ^)T/](x, 9) 

[a - a] (x, 9, x\ 9')Tf{x', 9')d^{x', 9') 



[& - 9, x\ 9')f{x' - r_(x', 9'), 9')d^{x', 9') 

\a-^]{x,9,x'^,9')f{x'^,9')d^i{x'^,9'). 

In the last equality above we change variables x'^ = x' — t-(x', 9')9' , and denote the 
distribution - -f]{x,9,x' - T-{x' ,9')9' ,9') := [a - a]{x,9,x' ,9'). 
On the other hand, from (!3T]) . we get 

[f{A - A)Tf]{x + r+(x, 9)9, 9) = [i« - + r+(x, 9)9, 9) 

_ + r4x, 9)9, 9, x'r, 9')fix'j„ 9')deix^, 9') 



R 

[a^ - a\x + r+(x, 9)9, 9, x'^, 9')f{x'ji, 9')d^^{x^, 9'). 



Therefore, in the sense of distributions 

(32) [a^ -a^]{x + T+{x,9)9,9,x' -T_{x',9')9',9') = [a - a]{x,9,x' ,9') 

Independently of the equality (!32|) above, due to the fact that a = a = in the 
Bji \ Q, by direct verification in the formula (ITTl) . we get 

(33) [A'^-A'^]{x'-t4x',9')9',9') = [A-A]{x',9'), {x',9')eT^. 
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Now ([32D and §^ applied to ^ yield for {x,e,x',e') eT+x r_: 

(34) 0^ - f]^]{x + T+{x,e)e,e,x' -T-{x',e')e',e') = - f3]{x,e,x',e'). 

The isometric relation ffTOl) now follows from fl32l) and flMl). □ 



3. The forward problem in three or higher dimensions 

In this section we recall the properties which define the albedo operator and its 
kernel's singular expansion when n > 3. 

Let T be the operator defined by the left hand side of ([1]) considered on Bji x S*""^, 
n > 3. From ([6]), the second and the third terms of T are bounded operators on 

[Br X S'"'"^) , while the first term is unbounded. We view T as a (closed) unbounded 
operator on L^{Bji x S*""^) with the domain 

D{T) = {ue L\Br X 5"~i); e-Vue L\Br X 5""^), n|r_ = 0}; 

see [12] . 

We work under either one of the following subcritical conditions that yield well- 
posedness for the boundary value problem ([T]) and ([3]): 



r{x,9) / k{x,9,9')d9' 



's 



'n — l 



< 1, 



(35) ess sup(2, gjgj^xS"-! 
or 

(36) a{x, 9)- I k{x, 9, 9')d9' > 0, a.e. n x 5""^ 

see, e.g., [3 [El [IS [21 [25] . 

Let 6{x}{x') represent the delta distribution with respect to the boundary measure 
dfj,{x') supported at x G dQ, and let 6{0}{9') represent the delta distribution with 
respect to d9 centered at ^ G S*""^. The following result is a recast of [T^ Theorem 
2.3] to the unit speed velocities. 

Proposition 3.1. Let n > 3. Assume that the direct problem is well-posed. Then 
the albedo operator A : L^{r^,d^) L^{r^,d^) is bounded and its Schwartz kernel 
a{x,9,x',9'), considered as a distribution on with {x',9') G parameters, is 
given by a = ai + a2 + as, where 

(37) a^{x, 9, x' , 9') =Mg_^e-/o^"'"" 

(38) a,{x,9,x\9') J<f\-^'\ P^^' '\-/;+'^'^*^' '^'a(.-s.,.)..^./^(.'+s.',.')rf. 

n[x) ■ 9 Jq 

X k{x' + t9', 9', 9)5{x>+te'+T+(x'+te',e)e}{x)dt 

(39) \n{x') ■ 9'\-^as GL°°(r!^; ^^(r^, d^). 
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4. Preliminary estimates 

This section extends a result from [7j for continuous coefficients to the class of 
essentially bounded coefficients as in This extension is crucial to the transporta- 
tion of the problem to a larger domain, if no boundary knowledge of the coefficients 
is available, and has the added benefit of simplifying the proof considerably. 

Although the presentation below concerns the unit speed velocity and measure- 
ments at the boundary, the results easily extend to the variable velocity and mea- 
surements in the rotating planes setting. The main novelty in this section is made 
possible by the following result on the approximation of the identity. 

Proposition 4.1. [T^ Theorem 8.15] Suppose (f G L^(]R") with J (f{x)dx = 1 and 
\^{x)\ < C(l + for some C,e > 0, and let <^t = t-Xx/t). If f e Lp(M"), 

for some 1 < p < oo, then f * (ft ^ f{x) with t — » 0^, for almost every x in the 
Lebesgue set of /- in particular, for almost every x G M", and at every x at which f 
is continuous. 

Corollary 4.2. There is a family of maps (pe^^^'^ G L^{T_,dC,), {xq,C,'q) G r_ and 
e > 0, such that, for any f G L°°{T^,d^) given. 



(40) lim /_ ^ey,,e',{x', 9')f{x', e')di{x\ 9') = f{x'„ 9',) 



whenever {xq,9q) is in the Lebesgue set of f. In particular, fHU]) holds for almost 
every {xq,9'q) G r_. 

Proof For (x'q, 9'q) G r„ and e > sufficiently small, let (x', 9') : U x V C M^"-^ ^ 
x'{U) X 9'{V) cdQx S"""^ be a coordinate chart with {x'q,9'q) G x'{U) x 9'{U) such 
that x'q = x'{0) and 9'^ = 9'{0). For {x',9') G r„, define 

D(u,v) 



(41) 4>e,.'„e',{x',9') 



\n(x') ■ 9' 



D{x',9') 



VMx'))Ve{vi9')), 



where (p{u) = l/(ci;„_i) for |n| < 1, ip{u) = for |n| > 1, and (Pei^u) = e~"'~^^ip{u/e). 
By ujn-i we denoted the volume of the (n — l)-dimensional unit ball. Then, for any 
e > 0, J ip^{u)du = 1 and 



<Pe,.'„9',{x',9')f{x',9')dax',0') = / y,,{u)^,{v)f{x'{u),9'{v))dudv. 

jR2n-2 

Apply the equality above to / = 1 to get \\(f>e,x'Q,e'J\L^{r+,d^) = 1- The conclusion 
follows from Proposition 14. 1[ □ 

We introduce the following notations. For {x',9') G r„, < t < r+(x', 6*') and 
9 G let 

(42) x+(t, x', 9', 9) := x' + t9' + r+(x' + t9' , 9)9, 

(43) Fit, x',9', 9) := e-/o"^''''^*'''''«(^'+*^'+^^'^)'^^e"^o 

denote the exiting point after one scattering at the point x'+t9' coming from direction 
9' into the direction 9, and, respectively, the total absorption along the broken line 
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due to this scattering. We note here that depends on the point x' + 16' where the 
scattering takes place and not exphcitly on t. This remark will be used later in fISB]) . 

Let [a, k), [a, k) be as in ([6]) extended by on 5^ \ Q, such that the forward 
problem is well-posed. All the operators bearing the tilde refer to (a, k) and are 
defined in a similar way with the ones for {a,k). For example A is the albedo 
operator corresponding to (a, k). Recall that n is the dimension of the space. 



Theorem 4.3. Let {a, k), {a, k) be as in 

following estimates hold: 
For n > 2, 



). For almost every {x'q,6'q) G r_ the 



(44) 



For n > 3, 

Jo Js"~^ 

(45) 



k - k\{xQ + te'o, e\ e)E{t, x'^, e'^, e)dedt < \\a - a\\+ 
+l|F-F||oo/ / k{xo + teie'o,e)d9dt. 



Proof. Let {xq,9q) e r_ be arbitrarily fixed and let (pe,x'f,,e'g G ^^0^-) be defined as 
in Corollary 14. 2[ To simplify the formulas, since {x'q,6'q) is fixed, in the following we 
drop this dependence from the notation 0^ = (pg^^^g'^. 

Let w4 = .4,1 + ^2 + ^3 be the decomposition of the albedo operator given by 



Aj{x,e) 



ai{x, e, x', e')f{x', e')d^i{x')de\ 2 = 1,2, 3, 



where a^, i = 1,2,3 are the Schwartz kernels in Proposition 13.11 and dfi{x') is the 
induced Lebesgue measure on the boundary dQ. 

Let (f) G L°°(T^) with ||0||oo < 1- Since ||0e||Li(r-) = I5 the mapping properties of 
the albedo operator imply that 



(46) 



^{x,9)[A- A](j)s{x,9)d^{x,9) 



< \\A-A\\. 



Next we evaluate each of the three terms in J^^ (j){x, 0)[A — Ajcpe^x, 9)dC,{x, 6) by 
using the decomposition in Proposition 13.11 and Fubini's theorem. 
The first term is evaluated using the formula ([3] 



(47) 

/l(0,£) 



<i){x,d)[Ai-Ai]<Pe{x,e)d^{x,d) 



ix' + T,ix\e'),e')d>Jx\e') 



■Jo 



a{x'+se',e')ds _ -/o'+^"' ' a{x'+se' fi')ds 



d^{x', 9' 
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Since the integrand above is in L°°(r_) by applying fj^Ol) . we get for almost every 



(4^ 



To evaluate the second term we use the notations (H2|) and (H3l) and the formula 



38]): 



X 



x', a;, ^)A;(x' + tO', 9\ 9) - x', 9', x, + t9' , 9', 9) 



dtd9 \. 



Apply again poj) for the continuous integrand above to obtain for almost every 



£^0 



(49) 



-F{t,x'„9'„x,9)k{x', + t9i9'„9) 

or /2(0) = /2,i(0) + ^2,2(0) with 
(50) 

/2,i(0) = 
(51) 

1^2,2(0)1 < 



0(x+(t, x'„ 9'„ 9), 9) [F{t, x'^, 9i X, 9)k{x', + t9'„ 9'„ 9) 
dtd9. 



0(x+(t, x'o, ^d, 9), 9)F{t, x'o, ^0, X, 9){k - k){x'^ + t^g, 9'^, 9)dtd9, 



'5"-i Jo 

Consider the third term 



\F - F\ (t, x[), ^(), X, ^)A;(x[) + 9)dtd9. 



(P{x,9)[A3-A3]Mx,9)d^{x,9) 
(t>e{x\9')di{x',9' 



By (139p . the map (x', 6»') t-^ 0(x,6')(a3 - a3)(x,6',x',6'')|n(x') ■ 6''| ^rf^(x, 61) is in 
L°°(r„), and then, by (140|) we get for almost every (xq,6'q) G r„ 

|a3 - a3|(a;,6',x'o,6'()) 



(52) /3(0)(x'o,e[,):=lim/3(0,£ 

£^0 



|^x,i 



Hx',)-9',\ 



-d^{x, I 
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The left hand side of fH6|) has three terms. We move the third term to the right 
hand side (with absolute values) and take the limit with e ^ to get 

(53) \IM + hmi^'o^O'o) < M + Um^'o^O'o), a.e. ix'„e',)eT., 

for any (p G L°°{T+) with ||0||oo = 1- 

We note that the negligible set on which the inequality above does not hold may 
depend on </>. In f l53|) . we shall choose two sequences of (p to conclude the estimates 
of the lemma. Since countable union of negligible sets is negligible, the inequality 
flS^ holds almost everywhere on r_, independently of the term in the sequence. This 
justifies the argument below for almost every {xq,9q) in r_. 

First we show the estimate (Hlj) . Let {(pm} be a sequence of maps such that 
|0m| < 1, 0m = 1 near (xq + t+{xq,9q)6q,9q) and with support inside {{x,9) G 
r+ : |x - x[) + r+(x'o, 9'o)9'f^\ + \9 - 9'q\ < 1/m}. Then (gHD gives 

(54) Ji(0„) = e-/o"+*''"'°^«K+«ei,e^,)d. _ ^- j;;+^^'o''>'o^ a{x',+se',,ei,)ds ^ 

independently of m. From pOl) we have \imm->oo l2{4>m) = since the support in 9 
shrinks to 9'q. From fl52|) we also have lim„_>oo -^sd^ml) = 0, since the support shrinks 
to one point. 

Next we prove the estimate fllSl) . Let {(pm,q} be a (double indexed) sequence defined 
by (j)m,g{x, 9) = xM¥>,ix, 9) where e L~(5"-i) with Xm = ^ for \9 - 9'^\ < 1/m 
and Xm = 1 for |6' — 1 > 1/m and ipq G L°°(r+) to be specified below. Regardless 
of the way we define ^Pq, the presence of Xm already yields /i(0m,g) = for all m, g 
since Xm{,x'Q + t+{x'q, 9'o)9'q, 9'^) = 0. 

Next we construct (pg. Note that we only need to have it defined for 9 ^ 9^ since 
the function x vanishes nearby 9^. 

Define first ipq{x,9) as follows. Using the notation (1421) . let 

Ti.'„e'„e := {x+(t, x'q, 9',, 9)edBn: 0<t< r+(x[„ 9',)} 

represent the set of the exit points of particles entering in the direction 9q at Xq which 
scattered once in the direction of 9. The set 11^' ai o lies in the intersection of the 
plane through Xq spanned by 9q, 9 and the boundary dB^. 
Now, for each 9 G 5""^ and x+(t, x'^, 9'q, 9) of (gj]) in the set 

(55) {x+ G dBn : distaB^(x+, n^^,e^,e) < 1/m}, 
we define 

(56) <^g(x+(t, x'^, 9'„ 9), 9) ■= sgn{k - ~k){x'^ + t9'^, 9',, 9). 

Outside the set in (155|) we define ipq{-,9) = 0. 

Since x~^(t,x'Q,9Q,9) in depends only on Xq + t^O'^o ^^'^ ^^e equation (ISBIl 
gives a well defined function G L°^{Br x 5*"^^). Note also that (13^ shrinks to a 
negligible set on c^S^^ with q oo. 

Now apply the estimate to and use h{(pm,q) = to get 

(57) |/2,i(0m,g)| (x[),^^;,) < \\A-A\\+Is{\(PmM<,^o) + |/2,2(0™,,)|(a;[„^[,). 
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Since the support of (j)m^g shrinks to a set of measure zero in r_|_ as g ^ cxd, we get 
for all m and almost every {xq,9'q) G r_, liiRg^oo h{\4>m,q\){xQ, 9q) = 0. Finally, from 
(l50l) we obtain for almost every (xq, 6q) G r_ that 

lim lim J2,i(0rrt,g) = / / F{t,x'Q,6'Q,x,6)\k - k\{x'Q + t6'Q,6'Q,6)dtd6, 

m— >oo n— >oo J gn—i J q 

while from flHTl) we have 

|/2,2(0,n,,)| < / / '\f - F\{t,x'o,e'o,x,e)~k{x'o + te'o,eie)dtde, 

Js"-^ Jo 

for all m, g. The estimate (145|) in the lemma follows. □ 



5. Stability modulo gauge transformations 



In this section we prove Theorem II .21 

We start with two pairs (a, k), (a, k) G t/s,p and let 

e := 

We shall find an intermediate pair (a', k') ~ (a, k) such that (fT3|) and (fT4l) hold for 
some constant C > dependent on S,p, and i?/?. 
Define first the "trial" gauge transformation: 

(58) ipix,e) := e-/o'"^""(«-'^)(-^W-, a.e. {x,9) eBrX 5""^ 
Then ip > 0, ip\^R = 1,9- V^(x, 9) G L°^{Br x 5""^) and 

(59) d{x,9) = a{x,9)-9-VJnip{x,9). 



Note that (/?|rfl is close, but not equal, to 1. We use the first estimate from Lemma |4]3] 
to decide how far from 1 can v^lr^ be. 

By OS]), we have for almost every (xq, ^q) G r_ 



< e. 



In each of the integrals of the left hand side above change variables t = t^{xq, 9^ 
and denote xq = Xq + r+(xo, 6*0)^0 S^^ 



— s 



(60) 



< e. 



When (xq,^o) covers almost everywhere we get that {xq,9'q) covers F;^ almost 
everywhere. 
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Now apply the Mean Value theorem to u e"" to get the lower bound 



(61) 



-Mo 



-uo 



T-ixo,e'o) 



(a — a)(xo — s9'q, 6'Q)ds 



ln(^(xo,^[,)| >e~2^^|ln(^(xo,e^, 



where Uq = Mo(xo, 6'g, a, a) is a value between the two integrals appearing at the 
exponent in the left hand side above, and ip is defined in (l58ll . 

From (160|) and (16T!) we get the following estimate for the "trial" gauge ip: 

(62) I In ip{x, 6) I < e^^^e, a.e. (x, 6) G T^. 

The "trial" gauge ip is not good enough since it does not equal 1 on T^. We alter 
it to some ip e L°°{Bn x S^) with 6 -Vlntp e L'^IBr x 5*""^) in such a way that 
'f'ldBn = 1- More precisely, for almost every {x,9) G Br x S"''^, we define (p{x,9) by 

r_(x,i 



(63) 



ln(^(x. 



In Lp{x, 



\nip{x + T+{x,e)e,e). 



Since < t_{x,9)/t{x,9) < 1 we get (p G L°°{Bji x S^). Following directly from its 
definition \yhp\qbr = 0- Indeed, for {x, 6) G F;^, we get that r_(x, 6) = Q and ip){x, 6) = 
1, whereas for (x, 6^) G F;^ we have that t_{x, 6) = r(x, 9) and x = x-|-r+(x, 9)9. Since 
both maps x i— t{x,9) and x i— ln(/?(x + t_^{x,9)9,9) are constant in the direction 
of 9 and since 6^ ■ V^r„(x,6') = 1, we get 6* ■ Vln(^(x,6') G L'^{Br x S"""^) and 

In (p{x + r+(x, I 



(64) 



9-V\n<p{x,9) = 9 ■V\nip{x,9) - 



t{x,9) 

Define now the pair (a', k') in the equivalence class of (a, k) by 
(65) a'(x, 9) := a(x, ^) - ^ ■ ln(^(x, 9), 



(66) 



A;'(x,^', 



(/?(x', 9') 



Now the albedo operator corresponding to (a', /c'), satisfies A' = A, and then 

(67) \\A' -A\\ = \\A-A\\=e. 

Next we compare the pairs (a', /c') with (a, /c) and show them to satisfy ( |T3l) and 

(ED- 

Using the definitions (jSH]), dnSD, the relation flM|) . and the estimate for ip on F;^ 
fl62|) . we have for almost every (x, 6*) G B^ x S**^"^: 

|a(x, 6*) — a'(x, 6^)1 = |[a — a](x, 6^) + [a — a'](x, 6^)1 

(68) = 1^ ■ In (^(x, ^) - e ■ In (^(x, e) | 



In (p{x + T+{x, 9)9, 9)\ ^ ^ e 



2RT. 



t{x,9) 



t{x,9Y 
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Since the coefficients are supported away from dB^ (by construction of Bji) such 
that flTUl) holds, following fIBSl) we obtain the estimate f[T^ in the form 



(69) 



\a- a Lo < £- 



Cr 



with cji from (fTOl) . 

Up to this point, all the arguments above also work for two dimensional domains. 
Next we prove the estimate (fT4|) . These arguments are specific to three or higher 
dimensions. Recall the formula ( H3l) adapted to a' 



F'(t,x',^', 



Jo 



and note that F = F' is a quantity preserved under the gauge transformation flB^ . 
This follows by direct calculation and the fact that = 1 on dB^. Then for almost 
all (x', 6') G r^, t e [0, r(x', ^')] and 6 G S'""^ we have the following lower bound 



(70) 



\F'{t,x',e',e)\ > e 



-4_RS 



Using the non-negativity of a and a' we estimate 
\[F-F']{t,x',9',9)\< 



J*a{x'+se',e')ds _ J^a'(x'+se' fi')dt 



+ 



< 



-/o 



a(a;'+t6l'+s6»,6»)ds 



< ee 



2_RS 



[a-a'](a;' + s^',^')^^s 
(is 



T+(x'+te'. 



+ 







T+(x'+i6l',6») 



[a-a'](a;' + t^' + se,e)ds 



ds 



(71) 



2_RS 



T{x' + sd\e') J, 

t(x' + t9': 



t 



t{x\ d> 



+ 



t{x' + tO' + s^, ^) 



The next to the last inequality uses (jHSD; the following equality uses the fact that both 
maps s ^ t{x' + s6', 9') and s h-> r(x' + t6'' + s9, 9) are constant in s, while the last 
inequality uses the fact that t < r(x', 9') and r+(x' + 16'' + 56*, 9) < t{x' + 19' + s9, 9). 
Therefore we proved that for almost all {x',9') G r_, t G [0,t{x' ,9')] and 9 G S*""^ 
we have 

(72) |[F-F'](t,x',^',^)| < 2£e2^^. 

Recall now the estimate (145|) with respect to the pairs (a', k') and (a, fc): 



(73) 



\k-k'\{x'o + t9Q,9',9)F'{t,Xo,9'Q,9)d9dt < ||^-^'|| + 



\F-F'\ 



k{xo + t9'Q,9'Q,9)d9dt. 
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Now use the lower bound for F' in ( ITOl) . the upper bound for ||F — F'||oo in (1721) and 
the hypothesis ||/i;||oo,i < P, to obtain 

(74) f [ \~k~ k'\{x'o + te'o, e', e)dedt < ee^^^ (1 + 2pe^^^) . 
Jo J s^-^ 

Finally, integrating the formula above in {x'q^O'q) G r_ with the measure (i^(xo, 6'o), 
we get 

(75) 11^- - k'W^ < ETiRe^'''^ (1 + 2pe2«^) . 
Now choose 

C = max{7ri?e^^^ (l + 2pe2^^) , e^^^/c^} 
with cr from ( fTOl) to finish the proof of Theorem 11.21 

6. Preliminaries for two dimensional domains 

This section introduces the framework for the problem in two dimensions. The 
results are mainly from [29] . 

As above, T denotes the operator defined by the left hand side of ([T]) in Br x 
with Br C M?. The coefficients are extended to be in Br \ Vt. By the regular- 
ity assumption ([7]), the second and the third terms of T are bounded operators in 
L°°{Br X S^). The first term is unbounded. We view T as a (closed) unbounded 
operator on L°°{Br x S^) with the domain 

D{T) = {ue L°°{Br X S^)] e-Vue L°^{Br X 5^), nirfi G L°^{T^)}. 

To simplify notation, for x ^ y, we denote by 

x-y 

x-y = arg(x - y) = r, 

\x - y\ 

the direction from y to x. Also, let 

(76) < E{y, x) = e- a{x-t{x-yy,^dt < ^ 

denote the attenuation along the segment in the direction from y to x. The attenu- 
ations corresponding to a' and a will be denoted by E' and E, respectively. 

The boundary value problem ([T]) and is equivalent to the operator equation 

(77) (/ - M)u = Jf., 
where, using (17^ . 

(78) Jf-{x,e) = E{x -T4x,e)e,x)f4x -T4x,e)e,e) 

(79) Kf{x,9)= [ k{x,e',e)f{x,9')d9', and 

POO 

(80) Mf{x,9)= E{x -te,x)Kf{x -te,e)dt. 

Jo 

Under the subcritical assumption 

(81) R\\k\\^ < 1/2, 
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the operator M : L°°{Br x S^) L°°{Br x S^) is contractive and ([77]) has a 
unique solution obtained by Neumann series. Moreover, for /_ G L°°{r^), we get 
u = {I — M)~^Jf^ G D(T) has a well defined trace in L°°{r^) given by 

(82) ^[u]{x,9):=[{I-Mr\jf^]\^H{x,9), {x,9)eT^; 

see [ni [22]. Therefore the albedo operator A : L°°(r^) ^ L°°(r^) is bounded and 
has the Schwartz kernel a{x,6,x' ,6') = </)|pfl(x, 6*; x', 6*'), where, for {x',6') G T^, the 
map (x, 6^) HH> 9; x', 6') is the fundamental solution of ([1]) subject to the boundary 
condition 

(83) 0|rH(-,-;a;',^') = |n(x') ■ ^^r'V}(-)V}(-)- 
More precisely, as shown in [291 Proposition 1], a = ao + «i + «2 with 

ao + «i + "2 = 700 + 701 + 702 := 0o|r« + ^0o|r« + ~ ^)~"^^^0o|r«) 
where 

rr+{x',e') 

(84) 00 = ^(x-r_(x,^)e,x)V}W / 6{x-x' -t9')dt, 

Jo 

(85) 01 = ^i^^Ml^i^ _ r4y, 9')9', y)E{y, y + r+(y, ^)^), 

(86) < 02 < rfi?||A;||^(l -ln|^' X ^1). 

The constant depends on R only, x is the characteristic function of Br^ and, for 
(x, 6*, x', 6*') G -Bi? X S*^ X r^, ?/ = ?/(x, ^, x', ^') is the point of intersection of the rays 
x' — > x' + oo9' and x ^ x — oo^. The kernel (3 in ([TB]) is then given by 

(87) /5(x, ^, x', ^0 = [701 + 702] (x, ^, x', ^0 , (x, 9, x', ^0 G x T^, 

where 7 is the trace operator on T^. 

The following estimate from |29] is needed later. 

Lemma 6.1. Let x be the characteristic function of Br, L{x',9') be the line through 
x' in the direction 9' and dl{y) be the Lebesgue measure on the line. Then 

(88) r xix - t9) [ I ^l^J M y)dt < C(l - In \9' x 9\), 
Jo Jl(x',0') \x-t9- y\ 

where C is a constant dependent on R only. 

7. Stability of the equivalence classes in two dimensions 

In this section we work under the hypotheses in Section [6] and prove Theorem 11.31 
Let (a, k), (a, k) G Vs,p be given with ||^ — ^||* = e. Define the pair (a', k') in the 
equivalence class of (a, k) by f[65|) and f[66|) as before, i.e. 

a'(x, 9) := a(x, 9) - 9 ■ V,. In (^(x, 9), k\x, 9', 9) := ^^^A;(x, 9', 9), 

where (p is given in (l64ll . 
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Then the corresponding albedo operator A' = A and, thus, 

(89) ll^'-^ll^ = ll^-^ll^ = e. 

In particular, 

(90) ||/3-/3'||oo<e. 



and ([60]) holds. 

Starting with fl60|l . the same arguments as the ones in the three dimensional do- 
mains, notably fl68|) . are valid in two dimensions to conclude the estimate fl69|) : 



(91) \\a-a'\\^<e{e^''^/cn)=:eC. 
In turn, flOTjl yields 



From the definition fl66l) and (p2l) we now get 

(93) ||A;'||oo<pe^^^^ 
Let 

(94) E^{y, 9', 6) := E{y - T.{y, 9')9')E{y, y + r+(y, 9)9). 

be the total attenuation along the broken path due to one scattering at y G Bfi, when 
coming from the direction 9' and scattering into the direction 9. The formula (1851) 
now reads 

(95) 0i(x, 9, x', 9') = xiy)[kE,]iy 0',9) ^ G i?^ x 5^ x r_, 

I (7 X (7 I 

where x{y) and y = y{x, 9, x', 9') are as described above in Section[6l We also consider 
E[ and defined similarly with the attenuation a' to replace a. 

The relation with the quantity in (|43]), F'{t, x', 9', 9) = E[{x' + t9' , 9', 9), allows us 
to use the estimates (1701) and (ITTl) to conclude 

(96) \E[{y,9',9)\ > e^^^^ {y,9',9) eBrxS^x S\ 

(97) ||^i-^;|U<25e2^^,. 

Now use (EED, (EZD, (ESD and (IHTl) in E[{k - k') = {E[ - Ei)k + {Eik - E[k'), 
(evaluated at (y, 9', 9) with y = y{x, 9, x', 9') as above,) to estimate 



e 



-4i?S 



^ - ^'1 < \E[ - Ei\k + l^iAi - E[k' 



< \E[ - Ei\k + IP - p'\ \9 X 9'\ + |70'2- 70211^ X 9'\ 
(98) < 2epe2^^ + £+ 1702 -70211^ x 9'\. 
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We estimate the last term using the identity 

\e' X ^||70'2 - 702| = \o' X e\-f{i - M')"^M'20[, - \e' x e\-f{i - m)-^m^^o 

(99) = \e' X e\-f{I - M')-^[M'^(f)'Q - M^4)l] 

+ \9' X ^|7(/ - M)-i[M' - M](J - M')~^M^0o 
To estimate the first term in the right hand side above we write 

[M'^(l)o - M^0o] = M'{M' - M)(j)'Q + (M' - M)M0[) + M2(0[) - 0o) 
and bound each of the terms as follows. From their definitions we have 



POO 

MM'(I)q = / E{x - to, x)dtx 
Jo 



x/ \ 77 \ E [x ,y)E {y,x -te)dl{y). 

L(x',e') \x — tu — y\ 



M' 



POO 

M'(I)'q = / E'{x-te,x)dtx 
Jo 



^ f k\x -te,x-te- y, e)k'{y, 6', x-te-y) , 

X / ■ ■ E [x ,y)E {y,x -te)dl{y). 

Jhix'fi') \x-tO - y\ 

Since (/ — M')~^ is bounded in with a norm dependent on the radius only, say 
C{R), by adding and subtracting one term and by using Lemma WT[ we estimate 

|7(/ - M'y\M - M')M'(f)'o\ <CiR)\\E - ^'||oo||A;||oo||A;'||oo(l - In \9' x 6) 

(100) + C{R)fk - A;'||oo||A;'||oo(l - In \e' x ^|). 
Now, from fl9T]) we get 

(101) \\E - E'Woo < 2R\\d - a'W < C{R, S, crJe, 

for some constant which only depends on R, S, cr. 

In what follows we keep the notation C{R, S, cr) for constants that may be differ- 
ent from equation to equation but they only depend on R, E, cr in an explicit, but 
inessential, way. 

Using the fact that < t(l - Int) < 1 for t E [0, 1], the bound < p and the 
bound in (193!) . we get from fllOOl) that 

(102) \e' X e\\-f{I - M'y\M - M')M'(f)'Q\ < C{R,E,cii){ep'^ + p\\k - k'\\^). 
By reversing the roles of M' and M, we get similarly 

(103) \e' X e\\-f{I - M'y\M' - M)M(t)'^\ < C{R,^,CR){ep'^ + p\\k - k'\\^). 
Similarly, from the definition of as above, we also get 

(104) 1^' X ^||7(/-M')"^M2(0[,-0o) < C(/?,S,cjf>p2. 
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The estimates ffT02|) . ffT03|) and ffTOij) imply 
(105) \e' X ^|7(/ - M'Y^[M'^<P'^ - M^^l] < C{R, S, cr){sp^ + p||^- - k'\ 



Next we estimate the second term of the right hand side of f l99l) . Since for any 
/ e L°°{Br X S^), we have 

[E'-E](x-t^,x) / k'{x-te,e',e)f{x-te)de'dt 



51 



E{x-te,x) / [k' -k]{x-te,e',e)f{x-te)de'dt 



\[M'-M]fix,9)\ < 
+ 

< |2i?||E' - ^llooll^'lloo + 2i?||A;' - A;||oo} 

< [C{R, S, ch)£p + 2i?||A;' - 

we get 

(106) \9' X ^1 |7(/ - M)-\M' - M)-'M^o\ < C{R, S, cr) [ep^ + \\k' - ^IUp'} • 
Since p < 1, by applying fllOSp and (11061) in flUU]) . we get 

(107) |7[02]-7[02]| \exe'\<C{R,T.,CR)[ep' + p\\k-k'\ 
Therefore the basic estimate (pH!) yields 

(108) II A;' - fc|U < C{R, E, cr)£ + C{R, S, c^)p|| A;' - ^lU- 
By choosing 

(109) p < -—^ -, 

C[R,Ts,cr) 

we get the final estimate 

mm 111.' Ill ^ C(fi.s.cB) , 

The constant C from Theorem 11.31 is the largest between the constant in ( I9T1) and 

dnoD. 

8. Concluding Remarks 

In the case of an anisotropic attenuating medium, the albedo operator determines 
the attenuation and scattering properties up to a gauge equivalence class. The set of 
gauge functions has a natural structure of a multiplicative group which acts transi- 
tively on the pairs of the coefficients. 

We showed that the gauge equivalent classes are stably determined by the albedo 
operator. We understand the distance between equivalent classes to be the infimum 
of the distances between the corresponding representatives. 

The proof uses essentially the fact that, without loss of generality, the problem 
can be transferred to a larger domain, and, consequently, the total travel time (with 
respect to the larger domain) of free moving particles in the interior domain stays 
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away from zero. The no loss of generality part is due to the extension of an estimate 
in [7] to essentially bounded coefficients. 

The fact that we get Lipschitz stability estimates in ( fT3l) . ( fT4|) instead of conditional 
Holder stability estimates as in [7] may seem strange. In fact, if we assume that a 
and a depend on x only (or on (x, |^|), if ^ belongs to an open velocity space), then 
(fT3|) implies 



in the L°° norm, compare with [3 Theorem 3.2]. Then, by using interpolation es- 
timates, and the stability of the X-ray transform, we can get a conditional Holder 
stability estimate for a — a, similar to the one in [TJ Theorem 3.4]. 



The third author thanks Jason Swanson for pointing out the approximation of the 
identity result in [T^ Theorem 8.15]. 
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